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ABSTRACT 


The work is devoted to the theoretical analysis of contact melting by the migrating heat source with 
an arbitrary shaped isothermal heating surface. After the substantiated simplification the governing cqua- 
tions arc transformed to the convenient for engineering calculations relationships. Analytical solutions 
are used for numerical prediction of optimal shape of the heating surface. Problem is investigeted for the 
constant and for temperature dependent physical properties of the melt. 


1. INTRODUCTION 


Melting of solids by contact with a heating surface takes place in numcro,^ natural and 
technological processes. These processes arc enumerated in the previous works [2, 4, 12-14, 22] devoted 
to contact melting problem and are divided into two groups. In one group the melting material lies on the 
heating surface and pressed against it by some external force (for instance, the force of the weight of the 
melting material). This situation arises when an unfixed solid melts in an enclosure [1, 16, 22] and in other 
contact melting devices used in industry [8]. Another group of applications involves a moving heat source 
melting its way through the surrounding solid. This situation arises in such fields as welding [21], geology 
[3], nuclear technology [9, 10] thermal drilling of rocks [4, 6, 18, 20] and glaciers [11, 17, 19]. Thermal dril- 
ling is commonly recognized now as the most effective method of boring glaciers. Boring rocks, sands and 
soil by thermopenctrators is a relatively new method in mining engineering. It has some advantages in 
comparison with traditional rotary drilling. The mo*t considerable advantage of thcrmodrilling is that 
three major facts of excavation (rock fracturing, debris removal and wall stabilization) arc accomplished 
in a single integrated operation. 

This work is devoted to the theoretical analysis of the contact melting process by the moving heating 
source with an arbitrary shaped isothermal heating surface. 

2. ANALYSIS 

2.1. The physical model and governing equations. 

Obviously every technological process where contact melting occurs has its own specific character. In 
particular case of thcrmodrilling. it is contact melting with a great specific load and heat energy, with ar- 
bitrary shaped heating surface. Thermopenctrators arc radially symmetric and in some cases ring-shaped, 
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^or toroidal, with a large central hole for forming and extracting the core sample [6], A schematic diagram] 
of the contact melting for the thcrmodrilling conditions is shown in Fig.l. Axisymmctric heater -1 is pen- 
etrating into the melting solid -2 with the velocity V under the effect of applied external force F. The 
thcrmopcnctrator is separated from the solid with a layer ofmclt -3, flowing along the thin channel be- 
tween the heating surface I h and solid-liquid interface £ m . It is assumed that the solid-liquid interface 
is a sharply defined surface and melting occurs precisely at temperature t m , melt flow is laminar and 
two-dimensional. Molten layer is assumed to be incompressible Newtonian liquid with a temperature de- 
pendent physical properties (except density). Experimental results [6, 15] indicated that the heat source ve- 
locity attains its quasi-steady, constant value V soon after initiation of melting. This fact justifies the next 
assumption of quasi-steady heat and mass transfer in the contact melting problem. 

According to the physical model and assumptions enumerated above the governing differential equa- 
tions of heat and mass transfer in the molten layer can be written as follows: 
dive = 0 


P L {v • V)e = G — p L Vp + divT 
c l p l (v' Vt L ) = div(;. L Vt L ) + d> 


where T is the deviator part of the tensor of internal stresses; v, p, t L arc the liquids velocity, pressure 
and temperature respectively; <D represents the dissipative terms in heal transfer equation; C L , p L , liq- 
uid properties defined in Nomenclature; the rest of the symbols arc standard. 

It is convenient for further analysis to use two systems of coordinates fixed to the healing surface: cy- 
lindrical coordinates (r, z) and local orthogonal boundary layer coordinates S and { are indicated in Fig. 1. 

Transforming (1) to non-dimensional form and using the similarity method in a preliminary analysis 
the main dimensionless parameters and numbers are generated [7]: 
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, =P L gd/W 


( 2 ) 


All the quantities here arc defined in the Nomenclature. Each of the dimensionless numbers (2) has 
an exact obvious physical meaning. In order to substantiate the simplification of the governing equations 
(1), the analysis of the values of these non-dimensional numbers for the concrete conditions of thermal 
drilling of ice and rock was carried out. Dimensionless parameter K h ~ 10~ 3 — 10" 1 physically represents 
the ratio of characteristic thickness of the molten layer and characteristic size of the heating surface d; cri- 
terion K g ~ 10 3 is the ratio of the characteristic mass force of the melt and external force; Reynolds 
number Re~ lO" 6 - lO' 4 ; Brinkman number Br~ 10" 5 - 10" 4 represents the viscous dissipation’ of heat in 
the molten layer, Pccklet number Pc~ 10—100; Stefan number Ste~ 1—10; K i( K c ~ 1. 

After ncglcctmg terms of 0(K b , K g , Re, Br) the governing nondimcnsional equations of heat and 
mass transfer in the molten layer will take the following form: 


Id. du 

— -(R Hu =0 

R • GtJ 

(3) 

as or\ of] 

(4) 
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where ij-jf P = p/w , S = s/d , H = h/K h d , /i = yi L //r“ , C = C L /C™ , /. 

_ v p vp t — t 

= A f /A f , R = r/d , u =-7 , u = , 0 = = — ,R = R(s)-cquation of generating line 

L L i Vp p Vp t — t 

» * h in 

T h of heating surface Y hi h = h(s) —thickness of the molten layer mersured along the internal normal to 
r h ; r 8 , v -longitude and transverse velocities in the molten layer; all the physical properties ofliquid C L , 

/jL, p L arc nondimensionalized by their values C^,/^,/x^at temperature t m ;rcfcrencc temperature ( 
t u — t ) is determined after nondimensionalization of Stefans condition, 

t -t = PcK (t -t ) ( 6 ) 

h m c m oc 


Here ^ is initial temperature of melting material. In the equation (3) v = 0 corresponds to the 
ring-shaped pcnctrator with a large central hole. In this case since the thickness of the liquid film is of 
0(K h ) it is possible to ignore the axially symmetric behaviour of heat and mass transfer and to consider (r, 
z) as the Cartesian coordinates; v= 1 corresponds to the continuous heating surface without hole. 

The boundary conditions in dimensionless form arc following 
At the heating surface Y h Of = 0) 

u, = u^ = 0; 0 — 0 h ; (7) 

0 L =(L—t )/(!.— t ); t is the unknown temperature on Y h 

b b m b in b 

At the solid— liquid interface Y. mfa = 1) 


u 


= 0 ; 


dR . 
ds ! 


0 = 0 


J _d0 

H dq * mi 


O dR 50 

=[ ^ +( lr )/Slc 1; Q=( “‘^r )l i„ 


( 8 ) 

(9) 


where 0 8 =( t 8 — t^ / (t m — t^), t g is a temperature of the solid material, n is an external relatively tomoltcn 
layer normal to Y m - 

For the pressure in the exit points of the molten layer s=s x and s = s 2 

P(Sj) = P(s 2 ) = 0 (10) 


Since when v = 1 it is only one exit point s= s : then in this case s l = 0 is the critical point where u, = 0 
dP / ds= 0. 

The assumption of quasi-stationary heat and mass transfer couscs the equality of external force F 
and the force of internal stresses in the molten layer. This condition with the defined accuracy of 0 (K h ) in 
non-dimensional form is [4] 

2 f R * 

1 ■ 2 J RPdR (11) 

r 2 - r ' J h, 


The function Q in Stefan condition (9) is the non-dimensional density of heat flux to the solid from 
surface Y m R/ s value depends oil the temperature distribution in the solid and is obtained from the solu- 
tion of the heat transfer problem which is the same as the problem of temperature distribution in the sur- 
rounding weldpool material [ 21 ]: 


L 


J 
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I SO c 0 i a 80 

V^ (R, ^ )=0; v=0 ’ ,; (12) 

0,l z =1; lim 0 =0; 0 1 \ =0; (R,Z) = (r,z)/ d; (13) 

0 f - is temperature distribution on surface £ f , formed after melting (Fig.l). Problem (12) (13) was 
solved in [4, 18] numerically by the boundary element method. 

2.2. Analytical solutions 

In [4] it was proved that boundary value problem (12), (13) admits an analytical solution as a func- 
tion of one independent variable when and only when the generating curve r m of £ is parabola. In 
parabolic coordinates a and t related to the coordinates R and Z by R = ter, Z= 0.5(a 2 — r 2 ) with boundary 
conditions on^ n ;i = T m ,fl 1 = 1, in infinity: t— oo, 0,— 0, cquasion (12) has the following solution 

Ei(-— t )/Ei(— — t ), v = 1 
2 2 m 

1 /pT (14) 

. crfc U T T) / c rre( / — T ), v = 0 


du x >0 


where erfc(x) = -2= [ cxp(-u 2 )du, Ei(-x)=f — — du x 

Vlt J , J , u 

According the formulae (9) and (14) the heat flux distribufion on is 


Pc c t 


f t + a 


— / i ~ a2 )’ v= 1 

'"fna. crfc(a), v = 0 


where a 2 = Pcr m / 2 Taking into account the fact that the distance between £ h and V m is the value 
of 0(K h ) we can rewrite (15) with the accuracy of 0(K t ) 

n D -« J dR fa 2 Ei( — a 2 ), v=l 

Q = Pc * c d 7 7 r- 06 ) 

ds *■ Vff acrfc(a), v = 0 

After simple transformation of equations (3) and (4) with invoked boundary conditions (7), (8) and 
(10) the velocities and pressure distribution in the molten Laver arc obtained 

R — R mtt — n 
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R . is a critical point which is determined by (19) and the boundary coudition P(s,) = 0 
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(20) 


when v = 1 it is supposed that R . = 0 

According the assumption 0 h = const the temperature distribution in the molten layer is sought as a 
function of one independent variable rj 

o = 0(n) 


As follows from intcrfacial condition (9) in this case 

— H = const, 
ds 

Last formulae (22) and (21) simplify equalities (18), (19) 

1 d R , , 

u = -=r -7-(P(r}) 

* D ds 
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— r 

h DdJ. 


r ,+, -r’ +1 


-ds 


(v+l)DJ, h 3 R' 

and heat transfer equation in the molten layer 

PcHtp(ry) dO _ d_ dO 
D dr] dr] ^drj 

Integration of this equation with the associated boundary conditions 


d 0 


0 u-°» E = 
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dR Ste 
ds 


reduce to the following relationship 

0(i;)=EHj 1 tcxp(^j ^drj)di7 

*1 " fl 

puting in (27) rj = 0 the temperature of the heating surface is determined 

„ — f‘l .PcH f 1 <pC , . 

8,-EKj jopf— J — d,Xl, 


(21) 

(22) 
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(28) 


In order to simplify further computations assume that r h is specified by the equation Z = A(R— R . ) 
In this case heat flux distribution Q is determined by the equality (16), where 

dR 1 


ds 


P(s) introduction into (11) yields 
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N /l + 4A 2 (R-R.) 2 
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(29) 


(Rj — R ,)(v + i) 

where R 1 = R.=0 if v = 1 

One of the most important characteristics of contact melting is the heat energy removal from the 
heating surface to the melt. Combining heat energy definition in non-dimensional form 


N = 2;rf J =(-;.^)| dR 
J R H v d 


with the equation (27) wc have 


L 
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r 


N = k(R j - R | )Ecxp(^p j —-dij) 


D J 


( 30 ) 




The quantity of heat energy calculated according (30) excesses the minimum heat power N 0 which is 
necessary to sustain the chosen melting velocity V. In non-dimensional form 

N o = n ( R 2 - R * K 1 + 1 /Stc) 

Here N 0 in comparison with N docs not contain the energy rate for heating melt and useless heat 
dissipation in the surrounding thcrmopcnctrator solid material. 

The main scope of present paper is to elucidate the influence of the heating surface shape upon the 
cITcctivcty of the contact melting process. Defining efficiency of the heating surface as a ratio /? = N a / N 
we'll have 


1 + 


/? = 


Ste 


exp( — PcH / 


D f^dn) 


(31) 


Equations (23), (24), (27)— (31) simulate heat and mass transfer processes in contact melting problem 
with the accuracy of 0(K h ). They arc convenient for prediction of contact melting process for materials 
with variable physical properties such as different kinds of rocks and sands. 

When the physical properties of melt arc constant (for example in the case of ice melting) equations 
(17), (23), (24), (27), (28), (30), (31). allows the considerable simplification, 

6(R’ +1 - R’ +1 ) 
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„ ( 1 + St^ ) , PcH, 

P = 5 exp( — ) 


(33) 

(34) 

(35) 

(36) 

(37) 


3. RESULTS 

Numerical prediction of u„ u„, P, H, 0, /? and other quantities of interest is carried out for ice and 
rock thcrmodrilling conditions. All the calculations of rock melting arc based on equations (23), (24), (27) 
— (31). Relatively complete description of basalt physical properties at high temperature is available in [5, 
6], Non-linear equation (27) is solved numerically by the iteration procedure. After this other quantities 
arc obtained automatically in a view of equations (23), (24), (28)— (31). As the initial estimate of iterative 
process solution (35) is chosen. When the ice boring process is investigated formulae (29), (32)— (37) arc 
used. The values of ice physical properties one can find for example in [4, 17, 19]. Effectiveness of the 
[heating surface is estimated by the value of parameter /?. Tt is shown in previous works [4, 17] that in com- 
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""parison with the other (cone-shaped, sphere-shaped, etc.) thcrmopcnctrators of the same power output 
the parabolic shaped pcnctrator attains the highest melting velocity, Therefore present paper is devoted to 
more detail analysis of contact melting with parabolic heating surface. The elongation of the surface is 
characterized by the value of shape parameter A. The results of numerically predicted efficiency as a func- 
tion of A for different conditions of ice and rock melting arc plotted in Fig. 2. Numerical results indicates 
that for slow melting when heat transfer in the molten layer is of minor significance and in the contrary 
the heat dissipation increases the flat heating surface (A < 1) is more cffccfivc. Vice versa for high speed 
melting the heat energy rate in the melt is dominating in comparison with the dissipation in surrounding 
solid material. So in this ease the elongate form of heating surface is preferable. According the calcula- 
tions presented in Fig. 2 there is the interval for melting velocities when the definition of the optimal shape 
is not trivial. For toroidal pcnctrator (v = 0)and ice melting conditions 20 < Pc< 55; for non— coring 
penetrator (v= 1) and rock melting conditions 40< Pc< 75. In order to find the maximum of P and the 
corresponding value of A, the derivative of p with respect to A is calculated. When the problem is 
non-linear and the physical properties of the melt depend on temperature the derivative is calculated 
numerically; when /i, c, / arc constants it is feasible to calculate P' A analytically. In a view of relationship 
(37) the equation P\ = 0 for computation of the optimal A can be written as follows: — PcH ' E — E' 

= 0 . 

This simple equation is solved by dividing segment in half method. 


L 


NOMENCLATURE 


A — shape parameter of the heating surface; 
a — thermal diffusivity; 
c — specific heat; 

d - characteristic size of heating device; 

F — external force; 

G - mass force; 

r, z - cylindrical coordinates defined in Fig.l; 
rj— internal radius of the heating device; 
r 2 — external radius of the heating surface; 
s lf s 2 — coordinates of the end points of generating curve of the heating surface; 
t — temperature; F-vclocity of the molten layer; 

v tt v n — longitude and transverse velocities in the molten layer; 

V - melting velocity; 

s, f- longitudinal and transverse local coordinates in the molten layer defined in Fig. 1; 

F 

W-spccific axial load from heating device side (W = — r — ) 

n ( r 2 ~ r ,) 

p— efficiency; 

T— generating curve of surface Y ; 

A- thermal conductivity; 

/i-dynamic viscosity coefficient; 
p- density; 

a, i— parabolic coordinates. 


g — acceleration; 
h — melt layer thickness; 
L — latent melting heat; 
P — pressure; 

Q — heat flux density; 
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Indices: 

L — liquid phase; h — heating surface; 

s — solid phase; m — melting surface; 

* — critical point; oo — value in infinite point. 

All the non-dimensional parameters, numbers and functions arc determined in the text: (2), (5), etc. 
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Fig. 1. Schematic representation of the contact melting proccss:l — heating device, 2~mclting solid, 
3— molten layer. 



Fig. 2. Ffficicncy as a function of shape parameter A. 

a) lee boring conditions; v = 0, R , = 3.2 (ring shaped pcnctrator) 

b) Rock boring conditions; v = 1 , R m = R t = 0 (non-coring pcnctrator) 





